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Abstract 



OO , We formulate the problem of efficient transport of a quantum particle trapped in a harmonic potential 

which can move with a bounded velocity, as a minimum-time problem on a linear system with bounded 
input. We completely solve the corresponding optimal control problem and obtain an interesting bang- 
bang solution. These results are expected to find applications in quantum information processing, where 
C^ ', quantum transport between the storage and processing units of a quantum computer is an essential step. 

They can also be extended to the efficient transport of Bose-Einstein condensates, where the ability to 
control them is crucial for their potential use as interferometric sensors. 



Index Terms 

Quantum control, quantum transport, linear systems, time-optimal control 

I. Introduction 

During the last decades, a wealth of analytical and numerical tools from control theory 
and optimization have been successfully employed to analyze and control the performance of 
quantum mechanical systems, advancing quantum technology in areas as diverse as physical 
chemistry, metrology, and quantum information processing [1]. Although measurement-based 
feedback control [|2l, ^ and the promising coherent feedback control [JH-^ have gained 
considerable attention, open-loop control has been proven quite effective since it does not require 
any quantum measurement, avoiding the associated problems. Controllability results for finite- 
and infinite-dimensional quantum mechanical systems have been obtained, clarifying the control 
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limits on these systems ||71-|l9l. Some analytical solutions for optimal control problems defined 
on low-dimensional systems have been derived, yielding novel pulse sequences with unexpected 
gains compared with those traditionally used lfT0l - ||23l . Numerical optimization methods, based 
on gradient algorithms or direct approaches, have also been used intensively to address more 
complex tasks and to minimize the effect of the ubiquitous experimental imperfections [|24ll - [|3T1l . 

At the core of modern quantum technology lies the problem of transfering trapped quantum 
particles between operational sites by moving the trapping potential. For example, most of the 
suggested architectures for the implementation of a quantum computer employ the transport of 
qubits from the storage to the processor unit and back, see ll32l . The transport should be fast 
and "faithfuU", i.e., the final quantum state should be equivalent to the initial one up to a global 
phase factor. Ideally, the absence of the vibrational excitations at the final site is required. The 
high-fidelity transport that satisfies this no-heating condition is characterized as frictionless. Note 
that frictionless quantum transport can be achieved by moving the trapping potential slowly in 
an adiabatic manner, where the system follows the instantaneous eigenvalues and eigenstates 
of the time-dependent Hamiltonian. The drawback of this method is the long necessary times 
which may render it impractical. A way to bypass this problem is to prepare the same final states 
and energies as with the adiabatic process at a given final time, without necessarily following 
the instantaneous eigenstates at each moment. The resulting final state is faithfuU while the 
intermediate states are not. 

This nonadiabatic regime, leading to shorter transport times [33], provides a privileged area 
for applying optimal control techniques. Numerical optimization methods have been used to 
calculate the optimal currents in a segmented Paul trap for fast transport of ions while sup- 
pressing vibrational heating [|34ll . [|35l . Fast quantum transport using optical tweezers, where the 
acceleration is altered in a bang-bang manner, has been demonstrated experimentally [36] and 
studied theoretically |[371 . For a moving harmonic potential, the limits of faithfuU transport with 
various types of imperfect controls have been evaluated [[38l . while an inverse engineering method 
using Lewis-Riesenfeld invariants has been employed to achieve efficient quantum transport in 
short times f3T\. 

In the present article, we study the problem of minimizing the time of frictionless quantum 
transport in the case of a harmonic trap moving with a bounded velocity. This is different from the 
case examined in ll36ll . where the acceleration rather than the velocity is bounded. A physical 
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system which can be modeled as a moving potential with bounded speed is the "magnetic 
conveyor belt" [l39l| . In this system, time-dependent currents in a lithographic conductor pattern 
create a moving chain of potential wells; atoms are transported in these wells while remaining 
confined in all three dimensions. The speed of displacement can be controlled by adjusting the 
frequencies of the modulating currents. In the next section we formulate the quantum transport 
with limited trap speed as a time-optimal control problem for a three-dimensional linear system 
with bounded input. Note that most of the examples presented in the literature are usually 
limited to two-dimensional systems, which allow the visualization of the optimal synthesis on the 
plane. The problem is completely solved in section [nil where an interesting bang-bang solution 
is obtained. The present study complements our previous work on minimum-time frictionless 
cooling of a quantum particle in a harmonic potential ll23l . 



II. Optimal Control Formulation of the Quantum Transport Problem 

The evolution of the wavefunction ?/'(x, t) of a particle in a one-dimensional parabolic trapping 
potential centered around the moving point s{t) is given by the Schrodinger equation [|40l 

dip 



%t 



+ ——{x-sY 



V', (1) 



2m dx'^ 2 

where m is the particle mass and h is Planck's constant; x is a scalar that varies on some compact 
interval and V' is a square-integrable function on that interval. We assume that the experimental 
setup is such that there are essentially no spatial restrictions due to geometrical constraints, for 
example the system is placed in the middle of a large enough vacuum chamber. When s(t) = 0, 
the above equation can be solved by separation of variables and the solution is 

oo 
lP{x,t) = Y,Cne-'^-'/'^n{x), (2) 



n=0 



where 



are the eigenvalues and 



E„ = ( n + ^ ) /j^, n = 0, 1, . . . (3) 



are the eigenfunctions of the corresponding time-independent equation 

+ -^—x^ ^„ = K^„. (5) 



2m dx"^ 
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Fig. 1 . Schematic representation of the frictionless atomic transport. The harmonic trapping potential is displaced by d, while 
the populations of all the oscillator levels n = 0, 1, 2, ... at the final time are equal to the ones at the initial time. The coefficients 
Cn should be independent of the spatial coordinate x. 



Here if„ in & is the Hermite polynomial of degree n. The coefficients c„ in ^ can be found 
from the initial condition 

/oo 
-oo 

Consider now the case where the trap is moving with a bounded velocity v{t) E [—V, V],V > 



s = v{t). 



(6) 



If s(0) = and s(T) = d, it corresponds to a displacement d of the system in the time interval 
[0,T], see Fig. [IJ For frictionless transport (no vibrational heating), the path s{t) should be 
chosen so that the populations of all the oscillator levels n = 0, 1, 2, . . . for t = T are equal to 
the ones at t = 0. In other words, if 



V'(x,0) = ^c„(0)^„( 



x] 



n=0 



and 



ilj{x,T) =J2cn{T)^n{x ~ d), 
n=0 

then frictionless transport is achieved when 

|c„(T)|2 = |c„(0)|2, n = 0,l,2, 



(7) 
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This is shown schematically in Fig. [TJ We emphasize that the coefficients c„ should be inde- 
pendent of the spatial coordinate x. Among all the paths s(t) that result in (|7]), we would like 
to find the one that achieves frictionless transport in minimum time T. In the following we 
provide a sufficient condition on s{t) for frictionless transport and we use it to formulate the 
corresponding time-optimal control problem. 

Proposition 1: If s{t), with s(0) =0 and s{T) = d, is such that the equation 

d + uj'^{a- s) = (8) 

has a solution a{t) with a(0) = 0, d(0) = and a(T) = d, d(T) = 0, then condition ^ for 
frictionless transport is satisfied. 

Proof: Without loss of generality we assume that the initial state is the eigenfunction 
corresponding to the n-th level ■ip{x,0) = '^n{x). We will show that when the hypotheses 
of Proposition [H hold then il){x,T) = e^'^"^'^^'^n{x — d), where 0n(T) is a global (independent 
of the spatial coordinate x) phase factor. This and the linearity of © imply that if ip(x,0) = 
Er=oCn(0)^„(x) then ^(x,T) = E^=oCn(0)e*'^"(^)^„(2; - d), thus condition © is satisfied. 
We follow Leach [.40,1 and consider the "ansatz" 

^(a;,t) = e*('^^+'^")*„(x-a), (9) 

where a{t) satisfies ([8]) and the accompanying boundary conditions, while 0„(t) is a function 
of time to be determined, with 0„(O) = 0. Observe that Q corresponds to a wavefunction 
centered around the moving point x = a{t). The choice of a phase linearly dependent on the 
spatial coordinate becomes physically transparent if we recall that the momentum operator is p = 
(h/i)d/dx [|4T]| . so the phase factor in ^ gives rise to an average momentum (p) = md. Note that 
because of the boundary conditions, we have i^i^x, 0) = \E'„(x) and ?/'(x, T) = e*'^"'^^^\E'„(x — d) 
for t = and t = T, respectively, therefore it suffices to show that ^ satisfies ©• Plugging ^ 
into (dJ, we obtain 



-{max + h4>n)'^ nix) 
where x = a; — «• Since 



md'^ h^ d"^ mu'^ 



+ — -(x-s)^ 



*n(x), (10) 



2 2m dx^ 
X- sf = X^ + '2{a- s)x + s'^ -a'^, (11) 
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(flOl ) becomes 



m[a + w^(a — s)]x + h(l)n 



— — + ^X^ + ^[a' + uj\s' - a")] \ ^M = 0. (12) 



2m dx^ ' 2 '^ ■ 2 
The coefficient of x is zero because of dS]). If we additionally use ([5]), then (fT2l) becomes 



m , 



h<Pn + En + -[a' + u\s' - a')] \ ^„(x) = 0. 



The following choice of 0„ 

0„(t) 



m 



-{Ent+— / [a^ + uj\s^ - a^)]c/t 



(13) 



(14) 



h r"" ■ 2 ,0 

assures that (fT3l) is satisfied, so dll) is a solution of ©. ■ 

We express now the problem of minimum-time frictionless transport using the language of 
optimal control. If we set (recall that V is the maximum trap velocity) 



u a uj , , V 

xi = —a, ^2 = y, a;3 = — s, u[t) = —, 



(15) 



and rescale time according to tnew = ^^^old' ^^ obtain the following linear system with bounded 
control, equivalent to ^ and ([8]), 



T 



where now x = (xi,X2,X3) and 



A 



X = Ax + u(t)b, 



1 

-1 1 





(16) 





1 



(17) 



The original transport problem is transformed to the following time-optimal control problem: 
Problem 1: Find the control u{t), \u\ < 1, which drives system (fT6l) from (0, 0, 0) to (7, 0, 7), 

'J = ud/V > 0, in minimum time. 

The boundary conditions on x are derived from those on a (see proposition [T]) and s. 



July 11, 2011 



DRAFT 



SUBMITTED TO ffiEE TRANSACTIONS ON AUTOMATIC CONTROL 7 

III. Time-Optimal Solution and Examples 

Before solving the optimal control problem, we establish the existence and uniqueness of the 
optimal solution using well known results for linear time-optimal processes. The following two 
theorems refer to the general linear system 

x = Ax + Bu{t), (18) 

where x eW, u e U CM."', A e M"^", and B = [bM . . . \bm] G M"^'" (hi G M"). 

Theorem 1 (Controllability of linear systems with bounded controls): Suppose that A is such 
that all its eigenvalues have real parts equal to zero. Let U be any control set that is a neigh- 
borhood of the origin in R™. Then the linear control system with controls in U is controllable 
whenever [JlZQ{A^hj,j = 1, . . . ,m} spans M" (theorem 6, chapter 5 in [|42l|'). 

Note that for the single-input case m = 1, like the system that we study in this article, the 
above theorem can be directly derived from the null controllability conditions. Recall that the 
sufficient conditions to be able to bring any initial state of a single-input linear system to zero 
(null controllability) are that the Kalman matrix has rank n, the control u = belongs to the 
interior of the control set, and the eigenvalues of matrix A satisfy Re (Aj) < [|43l . The full 
controllability requires additionally the null controllability for the system with matrix ~A, i.e.. 
Re (Aj) > 0, so that the original system with matrix A can be driven from zero to any final state. 
The requirements of theorem [H and especially that for Re (A,) = 0, are now obvious. 

Definition 1 (General position condition): Let the control set f/ be a convex, closed, and 
bounded polyhedron in M™. The matrices A, B, and the set U satisfy the general position 
condition if for every vector w, which has the direction of one of the edges of U, the vector 
Bw has the property that it does not belong to any proper subspace of R" which is invariant 
under the operator A; i.e., the vectors Bw, ABw, . . . , A^^^Bw are linearly independent. 

Theorem 2 (Existence and uniqueness for linear time-optimal processes): Let the control set 
f/ be a convex, closed, and bounded polyhedron in R'" satisfying, along with matrices A and 
B, the general position condition. If there exists at least one control which transfers the state of 
the system between two points, there also exists a unique optimal control that accomplishes the 
same transfer (theorems 13 and 11 in (111). 

Proposition 2 (Existence and uniqueness of the solution for problemU): Problem[T]has a unique 
optimal solution. 
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Proof: Matrix A in (fTTj ) has eigenvalues ±i, with zero real parts, while span{6, 
Ah, A%} = spanjea, 62, Ci} = M'^, where Cj are the obviously defined unit vectors. Additionally, 
the control set U = [—1, 1] contains the origin. From theorem \T\ we deduce that system (fT6l ) 
with u E U is controllable. So, there exists at least one control which drives the system from the 
initial to the final point. The general position condition is equivalent to the linear independence 
of vectors b, Ab, A%, which is true. From theorem |2] we conclude that there exists a unique 
optimal control that accomplishes this transfer. ■ 

Having established the existence and uniqueness of a solution, we move to solve problem [T] 
For a constant Aq and a row vector A G {E?)* the control Hamiltonian for the single-input linear 
system (fT6l) is defined as 

H = H{Xo, A, X, u) = Xo + X{Ax + ub). 



Pontryagin's Maximum Principle 11441 provides the following necessary conditions for optimality: 

Theorem 3 (Maximum principle for linear time-optimal processes): Let (x*(t), M*(t)) be a time- 
optimal controlled trajectory that transfers the initial condition a;(0) = xq of system (fT6l) into 
the terminal state x{T) = xt- Then it is a necessary condition for optimality that there exists a 
constant Aq < and nonzero, absolutely continuous row vector function A(t) such that: 

1) A satisfies the so-called adjoint equation 

A = -—- = -XA. 
ox 

2) For < t < T the function u 1— )■ H{Xq, A(t), a;*(t), m) attains its maximum over the control 
set U at u = U:t{t). 

3) H{Xo,Xit),x,it),u,it))=0. 

Definition 2: A control u : [0, T] — )• [—1, 1] is said to be a bang control if u{t) = +1 on [0, T] 
or u{t) = —1 on [0,T]. A finite concatenation of bang controls is called a bang-bang control. 

Proposition 3: For problem \T\ extremal controls are bang or bang-bang. The latter controls 
are 27r-periodic. 

Proof: For system (fT6l) with coefficients given by (fTTj ) we have 

H{Xo, A, X, m) = Ao + A1X2 + A2(x3 - xi) + A3M, (19) 
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and thus 

Ai = A2, (20) 

A2 = -Ai, (21) 

A3 = -A2. (22) 

Observe that H is a linear function of the bounded control variable u. The coefficient of u 
in iJ is $ = A3, the so-called switching function. According to the maximum principle, point 
2 above, the optimal control is given by m = sign$, if $ 7^ 0. From (|20l) and (|22|) we obtain 
Ai + A3 = c, a constant, so $ = A3 = c — Ai. Also, from (|20l) and (|2TI) we get Ai + Ai = 
(harmonic oscillator), and then Ai(t) = Asm(t + 9), where A and 9 are constants. Thus 

^{t) = c-Asm{t + 9). (23) 

The constants A and c cannot be simultaneously equal to zero since A = implies Ai = A2 = 
and c = Ai = implies A3 = 0, in contradiction with maximum principle which requires 
A = (Ai,A2,A3) 7^ 0. Thus the extremal controls are obviously bang or bang-bang, with the 
latter being 27r -periodic. ■ 

There is a simple way to visualize the extremal trajectories in two dimensions. It is based on 
the observation that the projections of these trajectories on the XiX2-plane are concatenations of 
trochoids. Recall that a trochoid is the locus of a point at some fixed distance from the center 
of a circle rolling on a fixed line. Indeed, if we set 

yi = xi- X3, ^2 = a;2 =F 1, 

for M = ±1, then we find 

yl + yl = constant, 

for each time interval where the control is constant, and 

2/1=1/2, (24) 

2/2 = -yi- (25) 

From the last equations we find that the angular velocity of the rolling circle is Wc = 1. The 
center of the circle is (X3, ±1), so the horizontal velocity is Vc = |x3| = |m| = 1 and the radius 
is Re = Vc/cOc = 1- The circle rolls without slipping on the line X2 = 0. In Fig. [2] we plot an 
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(a) First part 



(b) Second part 
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2 

(d) Total trajectory 



Fig. 2. Projection of the optimal trajectory on the 2;ia;2-plane for 7 = tt. The circles with center (xs, ±1) that generate this 
projection by rolling on the line 3:2 = with velocity ±3 = it = ±1 are also shown (a) The first part OA is a cycloid (b) The 
second part AB is a prolate trochoid, since A lies outside the rolling disc (c) The last part BC is symmetric to the first one 
(d) Total trajectory. 



extremal trajectory with two switchings, along with the rolling circles that generate it. The first 
part of the trajectory, OA'm Fig. |2(a)[ is a cycloid, since y^ + ?/| = 1 = i?^. The circle generating 
this part has center (0:3, 1) and rolls to the right since X3 = u = 1 > 0. When the control switches 
to u = —1, the center of the generating circle becomes (0:3, —1) and it moves to the left. Fig. 
|2(b)[ since now X3 = m = — 1 < 0. The corresponding trajectory part AB is a prolate trochoid, 
since the moving point lies outside the rolling disc. After the second switching, the center of 
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the generating circle becomes again (0:3, 1) and it rolls to the right, Fig. |2(c)[ generating again a 
cycloid BC. The total trajectory OABC is shown in Fig. |2(d)[ The trochoids which compose it 
are synchronized such that the center of the circle and the point (xi,a;2) arrive simultaneously 
at the points (7, 1) and (7,0), respectively, so the final point in R'^ is (7,0,7). Note that there 
is a symmetry between the initial and the final part. As we shall see later, this observation is 
the key for the optimal solution. 

In the next proposition, we use the geometric intuition developed above to calculate the system 
evolution under an extremal input. 

Lemma 1 (Main technical point): Let = to < ^1 < ^2 < ■ ■ ■ < ^n = 2^- The alternating 
control input 

u{t) = {-iy-\ tj^^<t<tj, j = l,...,n (26) 

drives system (fT6l) from the origin x(0) = to the point x(t„) with coordinates 



n-l 



Xi 



"i{tn) - x^{tn) = - sin(t„) + 2 ^(-1)-' ^ sin(t„ - tj) 

n-l 
X2{tn) - i-ir-' = - COs(tJ + 2 ^(-1)^-^ COs(t„ - t,) 



i=i 



(27) 
(28) 
(29) 



The control —u{t) drives the system to the symmetric point — x(t„). 

Proof: Note first that since ±3 = u, (|29l) is obvious for the input (l26l) . In order to prove ([27 
and (l28l) . we use the attached to the rolling circles "moving" coordinates y = (^1,^2)^5 where 



yi = xi - X3, 



y2 = X2- (-1 



li-i 



(30) 
(31) 



for tj_i < t <tj,j = 1, . . . ,n. Observe that in each time interval, yi and 7/2 satisfy the equations 
((24)) and (l25l) of the harmonic oscillator, so 



y{t-)=R{t,-t,^Mtti), 



where R{r) is the rotation matrix 



R{r) 



cos r sm r 



— sm T cos T 



(32) 



(33) 
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When the control switches, there is a discontinuity in 1/2 



j = l,...,n-l 



(34) 



expressing the change in the center of the generating circle, (xs, ±1) — )■ (xs, =f1). Note that 





2/(0^ 



from the definition (|30l ), ([3T|) of y and the initial condition a;(0) = 0. 
For n = 1 we find from (1321), (l33l) and (l35l) 



y(tr) = i?(ti-to)y(t+) = i?(ii)y(0+) 



— sin(ti) 

— cos(ti) 
thus (1271) and (|28l ) hold. Now suppose that they hold for n even (odd), so 

- sin(t„) + 2 Ei=i'(-1)'"' sin(tn - t,) 

cos(t„) + 2 E?=i (-1)'"' cos(t„ - t,) 



y{t'r 



But 



2/(^,t) = l/(i;) + (-1) 



n-l 




2 



y{K)T 




2 



(35) 



(36) 



(37) 



where the minus (plus) sign in (|J7l ) corresponds to n even (odd), and 

i/(t;+i) = i?(t„+i-t„)y(^n)- (38) 

Using dill), (HIl) and ([33]) in dMl) we find 

2/1 (^n+l) 

= - sin(t„+i - t„) cos(t„) - cos(t„+i - t„) sin(t„) + 
2{=Fsin(t„+i - tn) + ^(-l)^"^[sin(t„+i - t„) cos(t„ - t^) + cos(tn+i - t„) sin(t„ - t^)]} 



j=i 



n-l 



- sin(t„+i) + 2[(-l)"-^ sin(t„+i - t„) + ^(-1)^'-^ sin(t„+i - t,)] 



j=i 



- sin(t„+i) + 2 J^(-l)^ ^ sin(t„+i - t^ 
j=i 
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and 

l/2(t^+l) 
= Sin(tn+i - tn) Sin{tn) - COs{tn+l - tn) COs(t„ 



n.-l 



2{=Fcos(t„+i - tn) + ^(-1)-' ^[cos(t„+i - to) cos(t„ - tj) - sin(t„+i - t„) sin(t„ - tj)]} 

n-l 
= - COs(Wi) + 2[(-l)"-^ COs(t„+i - tn) + ^{-ly-^ COs(t„+i - t,)] 

i=i 

n 

= - cos(t„+i) + 2 ^(-1)^"^ cos(t„+i - tj). 

The induction step has been proved. To prove the last statement in the lemma we use the variation 
of constants formula ||45l for the linear system (fT6l) . which for a;(0) = gives 

x{t) = [ u{a)e^^'-''^bda. 







Obviously, the control —u{t) drives the system to the symmetric point — a;(t„). ■ 

Theorem 4 (Optimal solution): For the final point (7,0,7), with 2(p — l)n < 7 < 2p7r,p = 
1, 2, . . ., problem [T] has a unique optimal solution with 2p switchings 

u{t) = i-l)'-\ tj-i<t<t,, j = l,...,2p+l, (39) 

where the constant control time intervals Tj = tj—tj^i are such that the initial and final intervals 
are equal ri = r2p+i = r and are given by the solution of the following transcendental equation 

2r + 2(p - l)vr - 7 _ 2 ^^^-1 ( sinr \ ^ ^^^^ 



2p — 1 \2p — COST 

while the intermediate intervals are 

2r + 2(p - l)7r - 7 
r2k = ^^3Y ^ k = l,...,p (41) 

and 

2p7r — 2r + 7 

T2k+i = 27r - T2k+2 = 7: z , fc = 1, . . . , p - 1. (42) 

2p — 1 

The total minimum transfer time is 

_ 4p[r + (p - l)7r] - 7 
t2p+i ^^^ . (43) 

For 7 = 2p7r the optimal control is u{t) = 1 and ^2^+1 = 2p7r. 
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Proof: We study first the bang-bang extremals. Consider an extremal control of the form (|26l ) 
with 2p switchings. From lemma [T] we have that the final state satisfies the terminal condition 

x{t2p+i) = (7,0,7)^ when 

- sin(t2p+i) + 2 5^(-l)^-' sin(t2p+i - t,) = 0, (44) 

- cos(t2p+i) + 2 ^(-1)^"' cos(t2p+i - t,) + 1 = 0, (45) 

2p+l 

5^(-l)^-i(t,-t,_i) = 7. (46) 

i=i 

If we multiply (|44l) by i = v^^ and add (l45l) we obtain 

2p 

_git2p+i ^ 2 ^(-l)i-ie^(*2''+i-*^) + 1 = 0. (47) 

We express this relation using the constant control time intervals Tk = tk — tfc-i- Due to the 
sinusoidal form with period 2n of the switching function (l23l) . for a bang-bang control it is 

< Tfc < 2n, k = 1, . . . ,2p + I, as well as r^ + r^+i = 2n for A; = 2, 3, . . . , 2p — 1 and p > 2. 
Also r2fc are equal for k = 1,2, . . . , p and p > 2, while r2fc+i are equal for fc = 1, 2, . . . , p — 1 
and p > 3. Using these relations, the times appearing in (|47l) can be expressed as follows 

2p+l 2p+l 

t2p+l = XI (^'^ " ^'^-l) = XI ^'^ = ^1 + ^('^ " ^)^ + ^2P + ^2p+l (48) 

k=l k=l 

and 

'ii,^ '4^ f (2p - J - l)7r + r2p + r2p+i, j odd 

t2p+i - ti = X ^^^ ~ ^'^-i) = X ^-^^ = 1 • ("^^^ 

fc>i fc>j [ (2p-j)7r + r2p+i, j even 

Using (|48]) and ([491) in ([47]) we obtain 

_g«(n+T2p+T2p+i) _|_ 2r,g«(f2p+f"2p+i) _ 2pe*'^^''+^ + 1 = 

which leads to 

e-^ = ^^^-^. (50) 

2p — e*"^! 

By taking the absolute value on both sides in the above equation we obtain 

cos Ti = cos r2p+i ^ Ti = r2p+i or n = 2n - T2p+i. 
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Under the second choice, (|501 ) takes the form e^^^p = i which has no solution in (0, 2tt). So 

n = r^p+i ^ r, < r < 27r. (51) 

Using dSB, (EQI) becomes 

e-^ = e^^^ </> = tan-^ f-^^^^^ y (52) 

\2p — cos T J 

where the range of tan^^ is taken to be (— 7r/2, 7r/2). For < r2p < 27r, (|52l) implies 

T-2p = 20, < r < TT (0 > 0) 

(53) 
r2p = 27r + 20, tt < r < 27r (0 < 0) 

By expressing (l46l) in terms of r^ and using (|5TI) and the other relations for these time intervals 
we obtain 

2p+l 2p+l 

5^ {-\f-\t^ - tfc-i) = 5^ (-1)'^- Vfc = 2r - pr2, + (p - l)(27r - r^,) = 7, 

fc=l k=\ 

SO 

_ 2r + 2(p-l).-o^ 
^ 2p-l ^ ' 

Using dill) and dUl), ([53]) becomes 

/,(r) = 0, < r < TT 

fp{r) - 27r = 0, TT < r < 27r 



(55) 



where 



It is 



„ , , 2r + 2(p — l)7r — 7 . / sin r , 
/p^ = ^^W - 2 tan-i . (56) 

2p — 1 \2P ~ COST/ 



f'M) - 2 ut^ - ./rri'A , > /;(°) - »■ o<-<2- (57) 



1 2p COS r — 1 

2p-l 4p2 -4pcosr + 1^ 

Note that the above derivative attains its minimum value when the second fraction in the 

parenthesis is maximized. This happens for cost = 1, which maximizes the numerator and 

minimizes the denominator. From the above inequality we conclude that fp{r) is monotonically 

increasing in the interval (0,27r). Since, additionally, 

f (r\\ 2(p-l)7r-7 2p7r - 7 
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we observe that for 2(p — l)7r < 7 < 2pTc,p = 1,2, . . ., the equation /p(r) = has a unique 
solution in (0, vr). On the other hand, 

jpv ; 2p- 1 ' •^''^ ^ 2p-l 

Note that /p(7r) - 27r < for p = 1, 2, . . . and 7 > 0, while /p(27r) - 27r < for p = 2, 3, . . . and 
7 > 0. Only for p = 1 and < 7 < 27r it is /i(27r) — 27r > and then equation /i(r) — 27r = 
has a solution in (7r,27r). Comparing this with the solution of /i(r) = for < 7 < 27r and 
using (|48l ). (|5T| ) and (|54l) . we find that in both cases the total time is 4r — 7 so the latter solution, 
which lies in (0, vr), corresponds to a shorter path. 

Now consider an extremal control of the form (|26l) with 2p — 1 switchings. Working as above 
we find 

2p-l 

_e^*2p + 2 ^ (_i)i-ie^(<2p-<.) -1 = 0, (58) 



where now 



2p 2p 

t2p = 5^(tfc - tk-l) = ^ Tfc = Ti + 2(p - 1)tT + T2p (59) 

k=l k=l 



and 



2£^ ^ I (2p-j-l)7r + r2p, j odd 

t2p - ij = 2^(4 - tk-i) = 2^ ^fc = < • (60) 

k>j k>j [ (2p - J - 2)7r + r2p-i + r2p, j even 



Using (1591) and dgOj), dSSj) becomes 

_gi(Ti+r2p) ^ 2pe*^2p _ 2(p _ I)ei(r2p-l+r2p) _ 1 = Q, 

which leads to 

e^"i + 2(p - l)e"^^-' + e"^"^" = 2p. (61) 

Observe that this equality holds only if the time intervals ri, r2p and T2p-i are integer multiples 
of 27r, which is not the case since they take values in the interval (0, 27r). 

Next we consider the extremal control —u{t), where u{t) is of the form (|26|) with 2p switchings 
and < Tfc < 27r, A; = 1, . . . , 2p + 1. It is not hard to check that (I50l)-(l53l) remain valid, but now 



2p+l 2p+l 

J2{-l)''-\h - h-i) = Y,{-lf-\k = 2r - pr2, + (p - l)(27r - r^,) = -7, 
fc=i fc=i 
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and hence 

2r + 2(p-l )7r + 7 

2p-l 



^2p = ^^^^^ -. (62) 



Using (I52l) and ^^, (l53]) becomes 

^p(r) = 0, < r < TT 

5fp(r) - 27r = 0, tt < r < 27r 



(63) 



where 



Since 



, , 2r + 2(p - l)7r + 7 _i / sinr , 
^P^ = 7f T " 2tan . (64) 

2p — 1 V^P ^ COST ' 



^p(^) = f'pir) > 0, < r < 27r, (65) 

(yfp(r) is monotonically increasing in the interval (0,27r). However, since 7 > 0, it is also 

rn\ 2(p-l)7r + 7 ^ , . 2p7r + 7 

qJO) = -^ > 0, qJ-n-) = — > 

ypK J 2p - 1 ' ^''^ ' 2p-l 

and the equation gp^r) = has no solution in (0,7r). On the other hand, if r is a solution of 
the equation 5'p(r) — 2n = 0, n < r < 2n, then 

2r + 2(p - l)7r + 7 (2p + l)r + 2(p - l)7r + 7 



T2p + T2p+l = — h T 

> 



2p - 1 2p - 1 

(2p + l)7r + 2(p - l)7r + 7 (4p - l)7r + 7 



2p - 1 2p - 1 

Up - l)7r Up - 2)7r 

> ^-^ '— > ^-^ '— = 27r, 

2p - 1 2p - 1 

and the resulting control does not correspond to an extremal since it should be T2p + T2p+i < 2tt. 

For an extremal control —u{t), where u{t) of the form (|26|) . with 2p — 1 switchings and 
< Tfc < 27r, k = 1, . . . , 2p, it is not hard to check that (|6T1) remains valid, so there is no 
extremal control sequence of this form. 

Finally we study the bang extremals. The constant control u{t) = 1 drives the system to 
the points (2p7r, 0, 2p7r) at t = 2p7r,p = 1,2, .. . It is the only extremal control that achieves 
this transfer, thus it is time optimal. Note that u{t) = —1 drives the system to the points 
(— 2p7r, 0, — 2p7r) at t = 2p7r, p = 1, 2, . . . The proof of the theorem is now complete. Relations 
((4T]) and (gll) are derived using ^5^, while ([43]) is easily obtained from (gS]) using ([SB and (|54l) . 
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(a) Four switchings 



(b) Six switchings 



Fig. 3. Projections of the optimal trajectories on the 2;ia;2-plane for (a) 7 = 2.47r and (b) 7 = 4.47r. Solid (dotted) line 
corresponds to it = 1 (it = — 1). 




(a) 




(b) 



Fig. 4. (a) Minimum time i2p+i as a function of 7 G [0, IOtt] (b) t2p+i = t2p+i — 2(p—l)n as a function of 7 — 7 — 2(p— l)7r, 
7 £ [0, 2tt], where dashed (solid) line corresponds to p = 1 (p — )■ 00). 



In Fig. |3] we plot the projections of the optimal trajectories on the xia;2-plane for 7 = 2Ati 



and 7 = AAtt. In Fig. |4(a)| we plot the minimum time to reach the final point (7,0,7) as a 
function of 7 G [0, lOvr]. It is tempting to think that the plot segments for 2(p — l)7r < 7 < 2pn 
are translations of the initial segment (0 < 7 < 2n) by 2(p — l)7r in both axes. But this is not 
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the case, as we shall immediately see. Let 

7 = 7-2(p-l)7r, (66) 

hp+i = t2p+i -2{p-l)n. (67) 



The transcendental equation (|4()1) can be expressed as 



T — 7/2 1 / sinr , ^^„^ 

" tan-^ = 0, (68) 



2p — 1 \2p — COST 

while fl3] ) gives 



In the limit p — ;■ oo, (1681 ) becomes 



r- sinr -7/2 = 0, (70) 



and 



hp+i ^ 2r. (71) 

In Fig. |4(b)| we plot hp+i as a function of 7 G [0, 27r] for p = 1 (dashed line) and for p — )■ 00 
(solid line). The segments in Fig. |4(a)| approach the limiting case as p increases. 

IV. Conclusion 

In this article, we formulated the problem of efficient transport of a quantum particle trapped 
in a harmonic potential moving with a bounded speed, as a minimum-time problem on a 
linear system with bounded input. We completely solved the corresponding optimal control 
problem, obtaining an interesting bang-bang solution. Similar approach can be followed for 
the problem of atom stopping or launching. Additional restrictions on the control, reflecting 
possible experimental limitations, can be incorporated in the current analysis. The complexity of 
the resulting optimal control problems, which may increase the difficulty of the analytical study, 
can be overcome by using a powerful state of the art numerical optimization method based on 
pseudospectral approximations [|26l . [|29l . [f30ll . 

The results obtained here are expected to find application in quantum information processing, 
where quantum transport between operational sites is an indispensable step. They can also be 
immediately extended to the efficient transport of Bose-Einstein condensates [46], where the 
ability to control the condensates is crucial for their potential use as interferometric sensors 

m,m. 
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